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I. 


Statement  of  Work 


This  research  program  is  directed  at  developing  new  models  for  a  more 
fundamental  understanding  of  magnetism  and  to  investigate  new  magnetic 
phenomena  such  as  the  recent  discovery  of  exchange  coupled  multilayers .  The 
research  focuses  on  the  use  of  advanced  thin  film  preparation  techniques  to 
prepare  high  quality  multilayered  and  single  thin  films.  Much  of  the  research 
utilizes  magneto- transport  techniques  to  probe  the  magnetic  properties  of  the 
magnetic  films.  The  films  are  prepared  by  both  ultra  high  vacuum  sputtering 
and  molecular  beam  epitaxy  techniques.  In  particular,  the  magnetic  phenomena 
being  investigated  include  both  a  fundamental  understanding  of  the  spin- orbit 
coupling,  the  above  mentioned  exchange  coupling,  and  a  variety  of  interfacial 
effects.  The  interfacial  effects  include  exchange  coupling  between 
ferromagnetic  and  antirerromagnetic  films ,  electrical  transport  in 
multilayered  magnetic  systems,  and  the  effects  of  interfacial  morphology  on 
magnetization  pinning.  The  spin- orbit  coupling  effects  are  very  important 
from  both  a  basic  science  view  and  for  magnetics  applications.  In  terms  of 
applications,  the  spin-orbit  coupling  determines  the  crystalline  anisotropy 
energies  and  various  other  magnetically  important  quantities. 


II.  Magnetic  Properties  of  Magnetic  Multilayers  and  Thin  Iron  Films 

As  a  mechanism  of  describing  this  research,  we  will  list  and  discuss  the 
publications  which  have  been  either  published,  submitted  for  publication,  or 
are  currently  in  draft  form.  The  first  three  listings  and  reprints  were 
included  in  the  forecast  report  submitted  in  June,  1990.  The  fourth  was 
submitted  for  publication  and  has  since  been  published.  These  items  are 
included  for  completeness  in  this  annual  report. 


A.  "A  First  Order  Magnetic  Field  Induced  Phase  Transition  in  Epitaxial 
Iron  Films  Studied  by  Magnetoresistance,"  K.T.  Riggs,  E.  Dan  Dahlberg,  and  G. 
Prinz,  Phys .  Rev.  B41 .  7088  (1990). 


This  study  focuses  on  the  rotation  of  the  saturation  magnetization  in  the 
plane  of  the  epitaxial  iron  films  and  utilizes  the  anisotropic 
magnetoresistance  to  follow  the  rotation  of  the  magnetization  in  the  presence 
of  magnetic  fields  applied  parallel  to  the  plane  of  the  films.  As  mentioned 
in  A.  above,  the  surface  anisotropy  energy  in  the  (110)  films  is  insufficient 
to  rotate  the  magnetization  out  of  the  plane  of  the  films.  With  the 
magnetization  pinned  in  the  plane  of  the  film  then  the  rotation  of  the 
magnetization  in  the  plane  of  the  films  as  a  function  of  the  magnitude  and 
direction  of  a  magnetic  field  applied  in  the  plane  of  the  film  can  be  modeled 
as  a  first  order  transition.  The  easy  way  to  understand  this  behavior  is  to 
consider  two  easy  axes  separated  by  a  hard  axis  in  the  plane  of  the  film.  If 
the  magnetization  is  required  to  rotate  from  one  easy  axis  to  the  other  by  the 
application  of  a  magnetic  field,  the  applied  magnetic  field  must  be  of 
sufficient  strength  that  the  magnetization  can  pass  by  the  hard  axis.  Once 
this  occurs  the  magnetization  then  abruptly  (in  a  first  order  sense)  makes  the 
transition  to  the  other  easy  axis.  In  the  analysis  of  this  behavior  the 
uniaxial  and  fourth  order  anisotropy  energies  of  the  epitaxial  films  can  be 
determined. 


B.  "Magnetic  Domains  in  Epitaxial  (100)  Fe  Thin  Films,"  Jeffrey  M. 
Florczak,  P.J.  Ryan,  J.N.  Kuznia,  A.M.  Wowchek,  P.I.  Cohen,  R.M.  White,  G.A. 
Prinz,  and  E.  Dan  Dahlberg,  Mat.  Res.  Soc.  Symp.  Proc.  151 .  213  (1989). 

This  work  shows  how  the  surface  morphology  of  the  semiconductor  substrate 
can  influence  the  magnetic  properties  of  the  epitaxial  magnetic  films.  It  was 
found  that  the  misfit  dislocations  which  penetrate  the  surface  of  the  InAs 
alloyed  substrates  can  dominate  the  magnetic  properties  of  the  Fe  films.  The 
magnetics  were  investigated  with  both  a  magneto-optic  magnetometer  technique 
we  developed  (see  below  in  other  research)  and  with  imaging  Kerr  microscopy. 
The  most  interesting  feature  of  this  research  is  the  result  that  by  correct 
preparation  of  the  substrate  the  magnetization  can  be  made  almost  isotropic  in 
the  plane  of  the  film  and  controlled  with  modest  applied  magnetic  fields  (20 


C.  Magnetic  Anisotropy  Constants  of  Epitaxial  (110)  Fe/GaAs  Films  From 
77K  to  293K  Studied  by  Magneto-resistance,"  Daniel  K.  Lottis,  G.A.  Prinz,  and 
E.  Dan  Dahlberg,  Mat.  Res.  Soc.  Stop .  Proc .  151 .  213  (1989). 

Because  the  iron  films  are  locked  to  the  substrate  the  question  arises  as 
to  the  temperature  dependence  of  various  magnetic  properties .  In  this  paper 
we  studied  the  temperature  dependence  of  the  anisotropy  energies  and  K  .  A 
energy  is  also  found  in  bulk  iron  whereas  the  energy  is  unique  to  the 
epitaxial  iron  films.  This  paper  determined  that  the  K^  energy  is  consistent 
with  a  uniaxial  strain  arising  from  the  growth  of  the  metal  film  on  the 
semiconductor  at  elevated  temperatures.  A  comparison  of  the  K^  energy  shows 
differences  from  that  of  bulk  iron  but  the  origin  of  the  difference  is 
uncertain. 

D.  "Detecting  Two  Magnetization  Components  by  the  Magneto-optical  Kerr 
Effect,"  Jeffrey  M.  Florczak  and  E.  Dan  Dahlberg,  published  in  the  Jour.  Appl . 
Phys .  67 .  7520  (1990),  (reprint  enclosed) 

A  novel  technique  for  detecting  two  orthogonal  in-plane  magnetization 
components  was  developed.  This  technique  utilizes  the  magneto-optical  Kerr 
effects  to  sense  the  two  components.  These  components  of  magnetization  are 
parallel  and  perpendicular  to  the  plane  of  incidence  of  the  light  beam.  The 
ability  to  sense  two  components,  individually  or  simultaneously,  is  a  result 
of  the  disparity  in  the  longitudinal  and  transverse  Kerr  effects.  Based  on 
the  Fresnel  reflection  coefficients  of  these  two  effects,  an  analysis  is 
presented  describing  this  dual  component  sensitivity.  The  physical  conditions 
are  given  for  simultaneous  and  individual  detection  of  the  two  in-plane 
magnetization  components.  To  substantiate  this  analysis,  magneto-optical 
measurements  are  made  on  single  crystal  Fe  films.  The  results  are  discussed 
in  the  context  of  dual  component  sensitivity.  This  procedure  is  useful  for 
determining  the  magnetization  process  of  thin  films  and  as  a  probe  to 
determine  in  plane  preferential  growth  in  polycrystalline  films . 


E.  "Exchange  effects  in  MBE  grown  iron  films,"  Y-J  Chen,  D.K.  Lottis,  and 
E.  Dan  Dahlberg,  J.N.  Kuznia,  A.M.  Wowchak,  and  P.I.  Cohen,  accepted  for 
publication  in  the  J.  Magn.  and  Magn.  Mat.  (preprint  enclosed) 

The  films  prepared  in  the  MBE  system  are  allowed  to  form  a  passivating 
oxide  on  the  free  surface.  We  discover. d  that  this  oxide  orders 
antiferromagnetically  at  temperatures  on  the  order  of  200K.  A  comparisci,  of 
the  ordering  temperature  with  the  known  oxides  of  iron  allowed  us  to  identify 
the  oxide  as  FeO,  which  is  the  main  result  of  this  publication.  The  ordering 
of  this  oxide  is  manifest  by  its  effects  on  the  magnetic  properties  of  the 
underlying  iron  films.  The  hysteresis  loops  of  the  iron  films  exhibit  both 
an  increase  in  the  coercivity  of  the  iron  films  and  a  shift  in  the  hysteresis 
loops  due  to  the  exchange  coupling  to  the  antif arromagnetir  oxide.  Our 
investigation  of  this  phenomena  focused  on  the  temperature  dependence  of  both 
of  these  quantities.  The  research  on  the  exchange  coupling  between  the 
surface  oxide  and  the  iron  film  is  potentially  very  exciting.  At  the  present 
time  a  need  for  a  detailed  understanding  of  exchange  coupling  is  necessary  for 
several  technologies  from  magneto-optic  recording  media  to  bias  elements  in 
thin  film  magnetometers.  For  this  reason,  it  is  anticipated  that  research  on 
this  phenomena  will  continue . 


F.  "Simultaneous  in-plane  Kerr  effects  in  Fe/GaAs  (11C)  th  1  films,"  J.M. 
Florczak,  E.  Dan  Dahlberg,  J.N.  Kuznia,  A.M.  Wowchak,  and  P.I.  Cohen,  accepted 
for  publication  in  the  J.  Magn.  and  Magn.  Mat.  (preprint  enclosed) 

In  the  case  of  the  (110)  grown  films  which  have  a  hard  or  <111>  direction 
separating  the  easy  from  intermediate  directions,  the  magnetization  process 
occurs  via  a  discontinuous  jump  in  the  direction.  The  work  on  the  use  of 
magneto-optics  to  simultaneously  measure  the  magnetization  of  two  orthogonal 
components  of  the  magnetization  in  a  thin  film  (see  D.  above)  was  used  to 
study  how  the  magnetization  in  the  epitaxial  films  evolves  in  the  presence  of 
an  applied  magnetic  field.  It  is  this  technique  which  provided  the  data  we 
used  to  model  the  effects  of  surface  morphology  on  the  magnetic  of  the 


epitaxial  films  described  earlier.  In  this  publication,  the  technique  allowed 
the  observation  of  the  discontinuous  jump  in  the  magnetization  when  realigning 
past  the  hard  or  (110)  direction.  In  the  previous  work  we  were  not  able  to 
monitor  the  magnetization  continuously.  This  indicates  the  great  utility  of 
the  MO  technique. 

G.  "A  model  system  for  slow  dynamics,"  D.K.  Lottis ,  R.M.  White,  and  E. 

Dan  Dahlberg  submitted  for  consideration  of  publication  to  Phys .  Rev.  Lett, 
(preprint  enclosed) 

The  studies  of  the  dynamics  of  the  magnetization  process  in  the  thin  iron 
films  indicated  a  need  to  understand  in  greater  detail  the  nucleation  and 
magnetization  process  in  magnetic  systems.  In  this  paper,  we  address  how 
interactions  in  a  magnetic  system  can  give  rise  to  slow  dynamics.  We  used  the 
dipolar  coupling  to  provide  the  interaction  in  a  planar  system  of  spins.  The 
remarkable  result  is  that  even  when  treated  in  the  mean  field  limit,  the 
system  responds  logarithmically  in  time.  We  were  also  able  to  show  that  the 
model  replicates  a  stretched  exponential  over  six  decades  in  time.  The  most 
remarkable  feature  of  the  model  is  that  it  also  predicts  much  of  the  behavior 
observed  in  the  decay  of  the  remanent  magnetization  observed  in  high 
temperature  superconductors.  In  particular  the  quasilogarithmic  decay  and  the 
nonmonotonic  temperature  dependence  of  the  decay  slope.  In  general,  this 
paper  indicates  that  interactions  and  not  disorder  may  play  a  very  important 
role  in  a  number  of  physical  systems  from  structural  glasses  to 
superconductors . 


H.  "Magnetization  Reversal  in  (100)  Fe  Thin  Films,"  Jeffery  M.  Florczak 
and  E.  Dan  Dahlberg,  to  be  submitted  for  publication  to  The  Physical  Review. 

This  publication  is  within  a  week  or  two  of  being  ready  for 
publication.  It  also  utilizes  the  technique  we  developed  using  magneto- optics 
to  provide  a  simultaneous  measure  of  the  magnetization  in  two  orthogonal 


directions  in  a  magnetic  film.  This  work  focuses  on  (100)  oriented  iron 
films.  These  films  differ  from  those  we  investigated  earlier  in  that  they  do 
not  possess  all  thi.ee  primary  three  crystallographic  directions  in  the  film 
planes  but  instead  only  have  the  <100>  and  the  <110>  directions.  The  lack  of 
the  <111>  direction  completely  alters  the  magnetization  process.  The  (100) 
films  appear  to  magnetize  in  a  much  more  uniform  manner.  The  data  taken  on 
these  films  was  modeled  with  a  uniform  rotation  of  the  magnetization 
direction.  In  the  modeling,  we  were  able  to  determine  the  anisotropy  energies 
and,  as  expected,  show  that  prior  to  a  jump  in  the  magnetization  direction 
that  the  magnetization  process  is  dominated  by  wall  formation. 


I,  Near  term  publications:  In  addition  to  the  research  outlined  above, 
we  are  currently  studying  superlattices  of  alternating  layers  of  magnetic  and 
nonmagnetic  metals.  This  research  appears  to  be  providing  insight  to  the 
transport  properties  of  more  complex  systems  than  has  been  previously 
understood.  Our  use  of  the  anisotropic  magnetoresistance  to  determine  the 
fraction  of  the  scattering  in  the  magnetic  system  has  provided  the  key  to 
resolve  the  scattering  in  the  different  layers.  It  is  expected  that  there 
will  be  an  increase  in  this  area  of  research  as  our  expertise  increases.  The 
initial  samples  investigated,  prepared  by  Professor  I.K.  Schuller  of  UCSD  and 
Argonne  National  Lab,  were  of  Co  and  Ag.  Since  the  completion  of  this  initial 
work  we  have  started  a  similar  project  using  Fe  and  Cu  for  the  metals.  These 
samples  also  were  prepared  by  Professor  Schuller.  We  are  currently  preparing 
samples  in  our  UHV  sputtering  system  which  was  purchased  with  funds  provided 
by  the  University  of  Minnesota  and  an  AFOSR  equipment  grant  (AF0SR-89-0138) . 

Another  research  area  is  the  study  of  the  spin-orbit  interaction  by  the 
temperature  dependence  of  the  anisotropic  magnetoresistance  and  anisotropy 
energies.  In  the  films  we  have  grown,  it  is  anticipated  that  the  largest 
effect  of  an  altered  iron  lattice  constant  will  be  in  those  properties  which 
depend  up  on  the  spin-orbit  coupling.  This  has  been  shown  already  in  iron 
films  grown  on  GaAs  substrates  which  compresses  the  iron  lattice.  This 
compression  alters  the  spin-orbit  coupling  which  determines  the  anisotropy 
energies  in  magnetic  systems.  In  these  films  the  small  alteration  of  the 


lattice  constant  is  sufficient  to  rotate  the  easy  axis  from  the  (100) 
direction  to  the  (110)  direction.  The  spin-orbit  coupling  is  the  fundanu  ntal 
mechanism  which  determines  the  anisotropic  roagnetoresis tance ,  the 
magnetoelastic  coupling,  the  magneto-optic  coupling,  and  provides  a 
fundamental  limit  for  the  coercivitv  via  the  anisotropy  energies.  Although 
this  one  interaction  is  of  fundamental  and  technical  importance  it  is  still 
not  well  understood  theoretically.  For  this  reason  we  are  systematically 
studying  many  of  the  material  properties  which  depend  upon  the  spin- orbit 
coupling  in  order  to  provide  the  necessary  information  to  develop  a  basic 
understanding  and  model  of  this  interaction.  This  work  was  primarily  being 
investigated  by  D.  Lottis  a  PhD  candidate.  He  has  written  his  thesis  and  is 
currently  on  a  Postdoctoral  appointment  with  Dr.  A.  Fert  in  Orsay,  France  It 
is  anticipated  that  his  thesis  will  be  published  a  either  a  single  major  work 
on  the  above  topics  or  if  some  require  further  efforts,  it  may  result  in 
several  smaller  Physical  Review  papers. 
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Detecting  two  magnetization  components  by  the  magneto-optical 
Kerr  effect 
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Minnesota  55455 

(Received  Id  January  1990;  accepted  for  publication  26  February  1990) 

A  novel  technique  for  detecting  two  orthogonal  in-plane  magnetization  components  is 
presented.  This  technique  utilizes  the  magneto-optical  Kerr  effects  to  sense  the  two 
components.  These  components  of  magnetization  are  parallel  and  perpendicular  to  the  plane 
of  incidence  of  the  light  beam.  The  ability  to  sense  two  components,  individually  or 
simultaneously,  is  a  result  of  the  disparity  in  the  longitudinal  and  transverse  Kerr  effects. 

Based  on  the  Fresnel  reflection  coefficients  of  these  two  effects,  an  analysis  is  presented 
describing  this  dual  component  sensitivity.  The  physical  conditions  are  given  for  simultaneous 
and  individual  detection  of  tne  two  in-plane  magnetization  components.  To  substantiate 
this  analysis,  magneto-optical  measurements  are  made  on  single-crystal  Fe  films.  The  results 
are  discussed  in  the  context  of  dual  component  sensitivity. 


I.  INTRODUCTION 

Conventional  methods  for  obtaining  magnetization 
loops,  e.g.,  vibrating  sample  magnetometer  (VSM), 
SQUID  susceptometer,  and  BH  loopers,  are  often  con¬ 
structed  to  detect  only  a  single  component  of  magnetiza¬ 
tion.  A  hysteresis  curve  from  these  single-component  tech¬ 
niques  is  suitable  for  determining  quantities  such  as  the 
squareness,  remanent  magnetization,  and  the  coercive  field. 
These  parameters  are  important  in  the  characterization  of 
magnetic  materials.  A  hysteresis  curve  from  these  meth¬ 
ods,  however,  provides  only  a  limited  amount  of  informa¬ 
tion  about  the  reversal  process.  Additional  information 
about  the  reversal  can  be  obtained  by  rotating  the  sample 
in  an  applied  field.1  A  more  fundamental  method  for  in¬ 
vestigating  the  magnetization  process  entails  measuring 
multiple  components  of  the  magnetization.  Previous  work 
along  these  lines  has  proved  its  utility,  but  the  instruments 
are  rather  complex  to  cors*ri,ct.2,1 

In  this  paper  we  present  a  novel  method  for  detecting 
two  in-plane  magnetization  components  in  ferromagnetic 
materials.  This  method  utilizes  the  magneto-optical  Kerr 
effect  to  sense  the  magnetization  as  a  function  of  the  ap¬ 
plied  field.  The  magneto-optical  technique  is  not  only  sim¬ 
ple  to  implement,  but  can  be  used  to  probe  the  magnetiza¬ 
tion  in  small  regions  of  the  sample.  The  dual  sensitivity  is 
based  on  the  physical  distinction  between  the  transverse 
and  longitudinal  Kerr  effects.  The  longitudinal  effect  is 
characterized  by  a  rotation  of  the  plane  of  polarization;  the 
amount  of  rotation  being  proportional  to  the  component  of 
magnetization  parallel  to  the  plane  of  incidence.  The  trans¬ 
verse  effect  involves  a  change  in  the  reflectivity  of  the  light 
polarized  parallel  to  the  plane  of  incidence,  not  a  rotation 
of  the  polarization.  This  change  of  reflectivity  for  the  trans¬ 
verse  effect  depends  upon  the  magnetization  component 
perpendicular  to  the  plane  of  incidence.  These  differences 
are  described  in  Sec.  II. 

The  difference  in  the  two  Kerr  effects  can  be  exploited 
by  passing  the  reflected  light  through  an  analyzer.  For  a 


specific  case,  rotating  the  analyzer  alters  the  dependence  of 
the  transmitted  light  on  the  two  magnetization  compo¬ 
nents.  Hysteresis  curves  presented  in  the  discussion  section 
illustrate  this  connection  between  the  analyzer  orientation 
and  the  dual  component  sensitivity.  For  incident  light  po¬ 
larized  parallel  to  the  plane  of  incidence,  when  the  analyzer 
is  near  extinction  (minimum  light  intensity),  the  intensity 
depends  on  the  component  of  magnetization  parallel  to  the 
plane  of  incidence.  Rotating  the  analyzer  90°  from  extinc¬ 
tion  maximizes  the  response  to  the  component  of  magne¬ 
tization  perpendicular  to  the  plane  of  incidence.  At  angles 
intermediate  between  these  *wn  extremes,  both  compo¬ 
nents  are  simultaneously  detected  to  varying  degrees.  It  is 
this  degree  of  control,  coupled  with  the  simplicity  of  the 
technique,  that  makes  this  method  attractive  in  studies  of 
magnetization  processes. 

The  dual  component  detection  can  occur  in  both  poly¬ 
crystalline  and  single-crystal  films.  The  data  we  present 
here  represent  what  is  possiMe  with  single-crvstal  films.  In 
polycrystalline  films,  certain  conditions  must  be  present 
before  two  components  can  be  detected.  A  laser  “spot” 
comparable  to  the  dimensions  of  a  crystallite  in  a  polycrys¬ 
talline  film  can  result  in  the  detection  of  dual  components. 
This  would  be  useful  in  an  investigation  of  interactions 
between  crystallites  during  the  magnetization  process. 

In  the  next  section,  the  Fresnel  reflection  coefficients 
derived  from  a  phenomenological  theory  of  the  magneto¬ 
optical  Kerr  effects  are  listed  and  examined.  This  is  fol¬ 
lowed  by  an  analysis  describing  the  conditions  for  simul¬ 
taneous  or  individual  detection  of  the  two  magnetization 
components.  A  comparison  of  the  theory  and  data  for 
single-crystal  films  is  then  presented  in  Sec.  IV.  In  Sec.  V. 
some  conclusions  are  discussed. 

II.  THEORY 

The  magneto-optical  Kerr  effects  are  categorized  ac¬ 
cording  to  the  geometry  of  the  magnetization  in  relation  to 
the  plane  of  incidence  and  the  film  plane.  The  two  Kerr 
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FIG  I  Geometries  for  the  transverse  anti  longitudinal  Kerr  effects.  The 
orientation  of  the  analyzer  and  polarizer,  t>,„  are  measured  relative  to  the 
plane  of  incidence  A  clockwise  rotation,  as  viewed  toward  (he  lijtht 
source,  is  defined  as  a  negative  angle  while  a  counterclockwise  rotation  is 
a  positive  angle.  The  unit  vectors,  s  and  p  ,  denote  directions  perpen¬ 
dicular  and  parallel  to  the  plane  of  incidence,  respectively.  The  polariza¬ 
tions  of  the  incident  and  reflected  electric  fields  are  also  depicted 


effects  involving  in-plane  magnetization  are  depicted  in 
Fig.  1  In  the  transverse  Kerr  effect,  the  magnetization  is 
perpendicular  to  the  plane  of  incidence.  In  the  longitudinal 
Kerr  effect,  the  magnetization  is  parallel  to  the  plane  of 
incidence.  The  analysis  of  the  dual  component  sensitivity  is 
based  on  reflection  coefficients  derived  from  a  phenomeno¬ 
logical  model  of  the  Kerr  effects.4  In  this  model,  the  Kerr 
effect  is  described  by  a  dielectric  tensor  with  off  diagonal 
elements  These  elements  of  the  dielectric  tensor  depend 
upon  the  Voigt  magneto-optical  parameter  Q.  Residing  in 
this  parameter  is  the  interaction  of  the  electromagnetic 
field  with  the  magnetic  electrons.  By  virtue  of  this,  Q  is 
proportional  to  the  magnetization  in  the  ferromagnet.  The 
reflection  coefficients  listed  here  for  the  transverse  and  lon¬ 
gitudinal  Kerr  effects  are  only  approximate.  Terms  higher 
than  order  Q  are  neglected.  Typically  Q  is  small,  but 
there  are  exceptions.'  In  addition,  the  film  thickness  is 
large  compared  to  the  skin  depth  of  the  light  so  t hat  the 
substrate-film  interface  can  be  neglected. 

Starting  with  the  transverse  Kerr  effect,  the  Frcnsel 
reflection  coefficients  connecting  the  incident  and  reflected 
electric  fields  of  given  polarizations  are 


, 

/»(} 

rj')( 

h .  sin  2D 

rpr 

1  "fS  - 

PI  i 

J  _  _ 

( rr  cos  D  1  )  -  sm 

(  I  ) 

r\  -  (fl  n[))/(fi  v  ,,/P).  (2) 


In  the  above,  the  quantity  r’„x.  for  example,  is  the  coefficient 
for  the  transverse  effect  relating  the  incident  \  wave  to  the 
reflected  p  wave  The  v  and  p  refer  to  light  that  is  polarized 
perpendicular  and  parallel  to  the  plane  of  incidence,  re¬ 


spectively.  The  angle  of  incidence  measured  from  the  sam¬ 
ple  normal  is  0.  n  is  the  index  of  refraction  of  the  material, 
K:  =  1,1  Q >s  the  off  diagonal  element  of  the  relative  per¬ 
mittivity  tensor,  fi  -  cos  ft,  and  /i’  [1  (sin'  t))/n! j!  ' 

From  these  Frensel  coefficients,  one  should  note  that  the 
reflected  light  does  not  undergo  a  rotation.  The  only 
magnetization-dependent  quantity,  in  this  case,  is  the  re¬ 
flection  coefficient  connecting  the  incident  and  reflected  p 
polarized  light.  Practically  speaking,  this  implies  that  hys¬ 
teresis  loops  can  be  taken  without  a  polarizer  or  analyzer 
when  the  transverse  effect  is  present 

The  Fresnel  reflection  coefficients  for  the  longitudinal 
Kerr  effect  are 

r*pp=(nli  —  /?’)/( nft  t  f¥ ) .  (4) 

r[~((3  -  nP  )/(/?  +  nf]').  (5) 

/  /  YPK'- 

V  --  -  r-p  -  n -jp  ( nji  f  (f )  ( ft  ;  ~np ,  • 

where  the  superscript  /  refers  to  longitudinal  and  y  sin  0 
The  other  terms  are  defined  earlier  in  connection  with  Eq. 
(1).  The  above  coefficients  describe  a  rotation  of  the  inci¬ 
dent  polarization  by  the  presence  of  the  off  diagonal  terms, 
and  r1^.  In  contrast  to  the  transverse  Kerr  effect,  the 
diagonal  terms  here  do  not  depend  upon  the  magnetiza¬ 
tion.  A  final  note  about  all  the  coefficients  in  the  limit  that 
£-0,  implying  that  the  material  is  no  longer  magnetic,  the 
distinction  between  the  two  geometries  vanishes,  and  the 
reflection  coefficients  reduce  to  those  of  normal  metallic 
reflection. 

Generalizing  to  the  case  of  two  Kerr  effects,  consider 
plane  polarized  light  incident  upon  a  ferromagnetic  sample 
having  both  longitudinal  and  transverse  Kerr  geometries 
(see  Fig.  1.)  As  an  example,  consider  the  film  to  be  in  a 
single  domain  state  with  the  in-plane  magnetization  at  an 
angle  of  45"  relative  to  the  plane  of  incidence.  This  mag¬ 
netization  vector  can  be  decomposed  into  components  par¬ 
allel  and  perpendicular  to  the  plane  of  incidence.  If  the  film 
is  exposed  to  an  applied  field,  the  magnetization  will  alter 
its  orientation  to  minimize  the  energy  of  the  film.  As  the 
magnetization  changes  so  will  the  Frensel  coefficients  de¬ 
scribing  the  reflection.  Experimentally,  the  reflected  p  wave 
will  alter  in  amplitude  due  to  the  component  of  magneti¬ 
zation  perpendicular  to  the  plane  of  incidence.  The  varia¬ 
tion  of  the  v-wavc  component  in  the  reflected  light  is  a 
result  of  the  magnetization  component  parallel  to  the  plane 
of  incidence.  By  passing  the  light  through  an  analyzer  and 
rotating  it  so  that  either  the  v  wave  or  p  wave  is  transmit¬ 
ted,  thereby  blocking  the  other  orthogonal  component,  al¬ 
lows  selective  detection  of  either  component. 

The  remaining  analysis  focuses  on  the  relation  between 
the  intensity  detected  by  a  photodiode  and  the  polarizer- 
analyzer  angles.  First,  we  denote  the  transmission  axis  of 
the  polarizer  as  0  and  the  analyzer  as  0ir  For  an  observer 
looking  toward  the  photon  source,  the  angles  tire  defined  to 
be  positive  (negative)  for  a  counterclockwise  (clockwise) 
rotation  from  the  plane  of  incidence  (see  Fig.  I  ). 
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The  electric  held  after  transmission  through  the  polar¬ 
izer  has  the  vector  form 

K  ,  cos  fJ„  p  *  /f,  sin  ft.  s.  (7) 

In  this  expression  p  denotes  the  unit  vector  parallel  to  the 
plane  of  incidence  and  s  denotes  the  unit  vector  perpen¬ 
dicular  to  the  plane  of  incidence.  This  light  is  incident 
upon  tiic  ferromagnetic  him.  l  he  interaction  of  the  light 
with  the  multiple  Kerr  effects,  and  the  subsequent  relief - 
non.  can  he  expressed  in  terms  of  a  general  scattering 
matrix  for  two  in-plane  magnetization  component's  this 
matrix  is 

,S'":i  m :  .S'  (  Q‘  m . )  •  mSlQ/m j ) .  (  S  ) 

where  rn.  \t  ' M  and  m  M/M  .  Here  M  is  the  com¬ 
ponent  of  magnetization  perpendicular  to  t lie  plane  of  in¬ 
cidence,  M  is  the  component  of  magnetization  parallel  to 
the  plane  of  incidence.  V/  is  the  saturation  magnetization, 
and  the  scattering  matrices  for  the  transverse  and  longitu¬ 
dinal  Kerr  effects  are  SlQ/m.)  and  S(Q/m:).  respec- 
tivelv  I  liese  scattering  matrices  for  the  transverse  and  lon¬ 
gitudinal  Kerr  effects  can  he  reduced  to  matrices 
describing  only  the  reflection  at  one  interlace.  For  the  case 
of  the  transverse  effect,  the  matrix  is 

S'[Q/m.)  \  r"-'  (l>) 

\r,r  r  I 

The  term  in  parentheses,  Q/m„  denotes  that  wherever  a  Q 
appears  in  the  Frensel  coefficients  of  the  transverse  Kerr 
effect.  Q  m.  should  replace  that  Q.  Since  the  magnetization 
is  restricted  to  he  in-plane,  m]  +  mj  -  I .  Using  this  result 
in  relation  (8),  the  electric  field  after  reflection  from  the 
ferromagnetic  sample  is 

E'r={m]r’rp  -+  m]r'pp)Et,  cos  0p  +  sin  0p. 

E\  =  m)r\t ,£„  cos  0p  f  r[sEn  sin  0p.  ( 10) 
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This  reflected  lighi  is  sent  to  the  analyzer  where  onlv  the 
component  parallel  to  the  analyzer  transmission  ois  will 
pass.  Figure  2  show  s  the  decomposition  of  tills  eleci  i  k  tie  Id 
into  components  parallel  and  perpendicular  to  t fie  1  rails 
mission  axis.  This  transmitted  component  is  subsequently 
sensed  by  a  photodiode.  In  terms  of  the  reflected  elect ri,. 
field,  the  transmitted  light  is 

E.  -  E'  sin  IE  *  /;'  cos  IE.  i  1  1  i 

Substituting  Fquations  (  10)  into  Fq.  (11)  and  simplifying 
the  result,  the  electric  field  parallel  to  the  analyzer  trans¬ 
mission  axis  is 

E,=  £1,|  (»>,r'pr  +  m}r'rr) cos  0p  cos  0„  +  m'r.  sin  (0. 

-  0„ )  +  r'„  sin  0p  sin  | .  (12) 

The  output  current  from  a  photodiode  sensing  this  light  is 
proportional  to  the  modulus  squared  of  Eq.  (12).  Thus, 
the  resulting  expression  for  the  detected  signal,  normalized 
to  the  incident  intensity,  is 


///„=[ \m]r!pp  *-  mj/pp\:  cos2  cos 2  6,  +  |  2  shrffy  -  0(l)  +  \  r[s  \ 1  sin2  9p  sin2  6U  +  [  (m]r‘pp  +  m]Epr)m]r'p, 

4-  c.c.  |  cos  0p  cos  ()„  sin (8p  -  0a)  +  ( (m]rlpp  +  m]r'pp)r^  +  c.c.  ]  sin  9p  sin  0a  cos  6p  cos  9a 


-f  j  E/n'jr1’  +  c.c.  |  sin  0p  sin  0a  sin(0„  —  0„)  |. 


(13) 


III.  SINGLE-  AND  DUAL-COMPONENT  DETECTION 

In  the  experimental  setup,  due  to  the  use  ,,f  a  polarized 
laser,  it  is  convenient  to  fix  the  incident  polarization  and 
vary  the  analyzer  angle.  In  light  of  this,  it  is  desirable  to 
find  simplifying  cases  of  Eq.  (13)  for  fixed  polarizer  angles 
II..  Three  cases  are  found  that  reduce  expression  (13)  For 
two  of  these,  only  a  single  magnetization  component  is 
detected.  This  occurs  when  <)r  90“  (v  polarized)  and 

also  when  II.,  <)„.  When  <>r  90".  the  intensity  is  sensi¬ 

tive  to  variations  of  the  component  parallel  to  the  plane  of 
incidence  When  0  Ur  the  intensity  depends  upon  the 


component  perpendicular  to  the  plane  of  incidence.  For  the 
situation  in  which  0p  —  0:l  ~  90°,  a  possibility  that  appar¬ 
ently  satisfies  both  of  the  above  cases,  the  intensity  is  in¬ 
dependent  of  the  magnetization  and  therefore  neither  is 
observed.  The  third  simplifying  case.  0  -  0°  (light  polar¬ 
ized  parallel  to  the  plane  of  incidence),  results  in  the  de¬ 
tection  of  two  magnetization  components.  Because  of  the 
simplicity  and  importance  of  (his  sensitivity  to  both  mag¬ 
netization  components,  the  remaining  analysis  concen- 
nates  on  this  arrangement 

Substituting  II  0°  in  F.q.  (  13).  the  expression  simpli¬ 
fies  to 


7522 


,J  Appl  Phys  ,  Vol  67.  No  1 2,  1  5  Juno  1 990 


J  M  Florczak  and  E  D  Dahltierq 


7522 


Sample 


Optical  Beam  and 

Chooser  Solitter  Pnlnriypr  Mnnnfl 


FIG.  3.  Block  diagram  of  the  hysteresis  loop  tracer.  The  system  is  en¬ 
closed  in  a  thermally  stabilized  enclosure  to  minimize  mode  sweeping  of 
the  2-mW  HeNe  laser. 

///0=  |  m]r!pp  +  m]r‘pp | 2  cos2  9a  +  |  2sin2  9a 

-  [(w/'ip  +  m2:r‘pP)ml\r% 

+  c.c.  ]  cos  9a  sin  9a.  (14) 

This  relation  predicts  that  for  9a~G°  the  intensity  is  spe¬ 
cifically  sensitive  to  changes  in  the  magnetization  compo¬ 
nent  perpendicular  to  the  plane  of  incidence.  This  is  easy  to 
see  by  noting  that  relations  (1)  and  (4)  combine  in 
mjr1^  +  m1,Jpp  t0  give  ?  sum  of  a  nonmagnetic  term  and  a 
term  proportional  to  QMt/My  For  9a  =  90°  the  intensity 
depends  primarily  on  the  component  parallel  to  the  plane 
of  incidence.  This  component  appears  as  QM/MS  and 
|  Q\ 2MiM,/M2s  in  the  sin  9a  cos  9a  term  and  as 
Q'-Mj/M;  in  the  sin2  9a  term.  The  terms  quadratic  in  Q 
are  small  in  comparison  to  the  linear  terms  since  |Q|  i: 
taken  to  be  small.  At  analyzer  angles  between  these  ex¬ 
tremes,  the  intensity  depends  upon  both  components.  The 
degree  of  sensitivity  to  each  component  can  be  controlled 
by  rotating  the  analyzer  to  the  required  angle.  An  impor¬ 
tant  prediction  of  £q.('4)  concerns  the  last  term  on  the 
right.  The  sine  function  is  odd  in  the  angle  hence,  the 
sin  0,,  cos  0U  term  can  change  sign  depending  if  the  ana¬ 
lyzer  is  at  9 „  or  +  9 a.  The  distinction  between  positive 

and  negative  is  only  real  for  angles  0°  <  j  9U\  <  90”.  The 
significance  of  the  sin  9U  cos  9a  term  can  be  deduced  from 
the  data  by  comparing  the  +  0a  and  --  9a  hysteresis  loops. 

(V.  EXPERIMENTAL  DETAILS 

Figure  3  is  a  block  diagram  of  the  hysteresis  loop 
tr;-  The  system  is  maintained  in  a  thermally  stabilized 
b  a  to  minimize  mode  sweeping  of  the  HeNe  laser. The 
light  is  split  into  two  beams,  one  is  directed  to  the  sample 
and  the  other  is  sent  to  a  reference  photodiode.  The  light 
incident  upon  the  sample  is  polarized  m  the  plane  of  inci¬ 
dence.  Op  —  0°,  and  the  angle  of  incidence  is  maintained  at 
60°  in  al1  the  data.  The  intensity  detected  by  the  lock-in 
amplifier  is  A/  F{M)  -  Here  If.Vf)  is  the  signal 
from  the  sample  and  Iu.f  i,  the  intensity  from  the  reference 
beam.  The  voltage  from  the  reference  photodiode  is  ad- 
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justed  via  a  voltage  divider  so  that  A /  is  zero  when  the 
sample  is  demagnetized.  This  procedure  is  used  to  mini¬ 
mize  the  magnetization-independent  intensity  (normal  me¬ 
tallic  reflection  from  the  ferromagnet)  when  the  analyzer  is 
adjusted  from  extinction. 

The  analyzer  and  polarizer  are  both  Glan-Thompson 
prism  polarizers  which  have  an  extinction  ratio  of 
1  X  10  “  6. 10  This  high  extinction  ratio  is  important  for  sep¬ 
arating  the  s  and  p  polarized  light.  The  plane  of  incidence 
is  determined  by  adjusting  the  polarizer  and  analyzer  for 
minimum  intensity  when  the  film  is  demagnetized.  Near 
extinction,  the  intensity  depends  quadratically  on  the  po¬ 
larizer  or  analyzer  angle  so  that  the  plane  of  incidence  can 
be  determined  fairly  accurately." 

The  film  is  mounted  on  a  rotatable  sample  stage;  this 
allows  the  applied  field  to  be  oriented  at  any  particular 
angle  relative  to  the  crystalline  axes  of  single-crystal  films. 
The  external  field  is  applied  in  the  plane  of  the  film,  and 
perpendicular  to  the  plane  of  incidence.  The  intensity  data 
can  be  stored  digitally  or  with  an  X-Y  recorder  for  later 
analysis.  The  hysteresis  curves  are  typically  acquired  in  less 
than  2  min  and  at  room  temperature. 

V.  RESULTS  AND  DISCUSSION 

To  support  the  above  analysis,  measurements  were 
made  on  single— crystal  Fe  films  epitaxially  grown  on 
( 100)  GaAs.  The  preparation  of  these  samples  have  been 
previously  discussed.12  The  epitaxial  growth  of  Fe  on  the 
(100)  plane  of  GaAs  constrains  the  Fe  to  form  a  (100) 
plane.  In  the  (100)  plane  of  Fe  there  are  two  magnetically 
easy  (100)  axes  and  two  magnetically  intermediate  (1 10) 
axes.  The  remaining  discussion  focuses  on  the  hysteresis 
curves  taken  with  the  applied  field  nearly  parallel  to  a 
magnetically  intermediate  (110)  axis.  It  should  be  empha¬ 
sized  that  all  the  data  are  taken  with  the  same  orientation 
of  the  applied  field,  only  the  analyzer  angle  is  varied. 

This  first  example  centers  on  the  measurement  of  the 
magnetization  component  parallel  to  the  applied  field.  The 
previous  analysis  indicates  that  for  (9„~0°  any  intensity 
change  is  due  to  the  component  of  magnetization  perpen¬ 
dicular  to  the  plane  of  incidence.  The  hysteresis  curve  for 
9j  =  0°  is  depicted  in  Fig.  4(a).  With  reference  to  Fig.  1.  in 
this  experimental  setup  the  component  of  magnetization 
perpendicular  to  the  plane  of  incidence  is  also  parallel  to 
the  applied  field.  As  the  field  is  increased  this  component 
grows  in  a  manner  indicative  of  coherent  rotation.  When 
the  field  surpasses  approximately  2A’1//V/S  =  550  Oe  (bulk 
Fe  constants),  the  magnetization  saturates  in  the  direction 
of  the  applied  field.  Hysteresis  curves  obtained  from  a 
SQUID  susceptometer  with  the  field  parallel  to  the  mag¬ 
netically  intermediate  axis  of  this  sample  are  similar  to  this 
curve.  Thus,  for  this  analyzer  orientation  this  method 
could  be  directly  compared  to  other  single-component 
measurement  techniques. 

The  second  example  corresponds  to  sensing  the  mag¬ 
netization  component  perpendicular  to  the  applied  field. 
Figure  4(bl  depicts  the  hysteresis  curve  for  6,-90°.  At 
this  analyzer  orientation  the  intensity  variations  are  a  re¬ 
sult  of  the  magnetization  component  parallel  to  the  plane 
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FIG.  4.  Hysteresis  curves  for  (a)  0„  =  O°  and  (b)  6„=  —  85"  when  the 
applied  field  is  along  a  (110)  magnetically  intermediate  axis  of  the  crystal. 
The  9a  =  0°  data  is  scaled  to  bring  the  +  M,  and  —  A/,  saturation  states 
to  arbitrary  intensities  of  +  1  and  -  1,  respectively.  The  scaling  of  (b) 
and  the  curves  in  Fig.  5.  are  relative  to  this  0„  =  O"  curve. 


of  incidence.  By  similar  reasoning  as  before,  the  magneti¬ 
zation  component  parallel  to  the  plane  of  incidence  is  also 
perpendicular  to  the  applied  field.  As  the  field  is  increased 
and  the  magnetization  gradually  increases  in  the  direction 
of  the  applied  field,  the  component  perpendicular  to  the 
applied  field  vanishes.  This  is  true  regardless  of  the  polarity 
of  the  applied  field.  This  explains  why  the  light  intensity  is 
so  similar  at  the  two  saturated  states.  At  low  fields  near  0 
Oe,  due  to  the  orientation  of  the  ( 100)  easy  axes,  the  mag¬ 
netization  is  angled  45°  from  the  plane  of  incidence  This 
maximizes  the  component  perpendicular  to  the  applied 
field  for  this  sample.  Hence,  the  intensity  values  have  a 
maximum  difference  near  0  Oe  between  the  positive  and 
negative  magnetization  directions.  The  reversal  field  for 
this  magnetization  component  is  rather  large  compared  to 
the  component  parallel  to  the  applied  field.  250  versus  20 
Oe.  This  disparity  is  a  result  of  the  particular  magnetiza¬ 
tion  process  occurring  in  the  film.1'  With  this  orientation 
of  the  analyzer,  d:I~  90°,  it  is  possible  to  determine  the 
changes  of  the  magnetization  component  parallel  to  the 
plane  of  incidence. 

In  this  last  example,  we  examine  the  magnetization 
loops  resulting  from  the  detection  of  two  magnetization 
components.  The  simultaneous  detection  of  two  magneti¬ 
zation  components  is  predicted  to  occur  at  I),  angles  be¬ 
tween  0°  and  90°.  Depleted  in  Fig  5  are  hysteresis  curves 
for  analyzer  angles  of  j-  60°  and  -  60'.  A  comparison  of 
the  positive  and  negative  analyzer  angles,  •  0 :  and  0a, 
determines  those  regions  w  here  the  .sin  fljeos  0  term  of 
relation  ( 14)  s  if  importance.  For  this  sample,  the  sign 
dependent  term  is  significant  for  fields  less  than  about  550 
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FIG.  5.  Hysteresis  curves  for  Ihe  same  applied  field  orientations  as  in  Fig. 
j  but  with  analyzer  angles  of  (a)  0U=  -  60”  and  (b)  0a  =  +  60°.  The 
intensity  variation  here  depends  upon  both  magnetization  components. 
The  first  transition  ccurs  at  a  field  of  20  Oe  while  the  second  transition 
occurs  at  a  field  of  250  Oe. 


Oe.  The  simultaneous  detection  of  both  in-plane  compo¬ 
nents  is  inferred  by  the  two  transitions  appearing  in  these 
data.  The  first  transition,  relative  to  0  Oe,  occurs  at  20  Oe 
and  the  second  at  250  Oe.  These  reversal  fields  are  similar 
to  the  reversal  fields  observed  in  Fig.  4  for  the  components 
parallel  and  perpendicular  to  the  applied  field,  respectively. 
The  overshoots  and  undershoots  in  these  curves  are  a  re¬ 
sult  of  the  particular  mixing  of  the  two  Kerr  effects  and  are 
further  evidence  of  the  dual-component  detection.1' 

We  now  turn  to  the  question  if  this  technique  can  be 
used  on  polycrystalline  films.  The  above  analysis  assumes 
the  laser  illuminates  a  single  crystal.  The  results  are  still 
valid  if  the  area  illuminated  by  the  laser  is  a  single  crystal¬ 
lite  of  a  polycrystalline  film  or  a  polycrystalline  film  that 
has  an  in-plane  preferred  growth  direction.  In  an  isotropic 
polycrystalline  film,  if  the  area  illuminated  contains  many 
crystallites,  a  superposition  of  the  individual  magnetization 
vectors  in  each  crystallite  results  in  a  magnetization  vector 
parallel  to  the  applied  field.  Hence,  for  a  large  laser  spot 
two-component  detection  would  not  be  expected.  This  was 
confirmed  by  magneto-optical  measurements  taken  of  poly¬ 
crystalline  Fe  films.  The  shape  of  the  hysteresis  curve  did 
not  depend  on  the  analyzer  angle  as  in  the  single-crystal 
films. 

For  the  case  of  a  polycrystalline  film  with  an  in-plane 
preferred  growth  direction,  many  of  the  crystallites  will,  to 
some  degree,  have  their  crystalline  axes  aligned.  A  super¬ 
position  of  the  magnetization  vectors  in  this  case  can  result 
in  a  net  magnetization  with  a  component  perpendicular  to 
the  applied  fielu  t  he  detection  of  two  components  in  this 
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situation  would  be  similar  to  that  of  single-crystal  films. 
This  is  where  these  measurements  could  serve  a  dual  pur¬ 
pose  by  detecting  any  in-plane  anisotropy  of  polvcrystal- 
line  films  due  to  growth  conditions  or  substrate  texture. 

VI.  CONCLUSIONS 

This  paper  presents  a  novel  technique,  based  on  the 
magneto-optical  Kerr  effects,  for  monitoring  two  compo¬ 
nents  of  the  magnetization  as  a  function  of  the  applied 
field.  For  the  particular  situation  when  the  incident  light  is 
polarized  in  the  scattering  plane,  a  rotation  of  the  analyzer 
to  90°  or  0°  detects  magnetization  components  parallel  or 
perpendicular  to  the  plane  of  incidence.  At  angles  interme¬ 
diate  between  these  extremes  both  components  are  simul¬ 
taneously  sensed  to  varying  degrees.  With  suitable  modifi¬ 
cations  of  the  beam  diameter  of  the  laser  or  anisotropic 
growth  conditions  for  films,  this  technique  can  be  used  on 
polycrystalline  films. 

The  in-plane  magnetization  components  can  also  be 
detected  simultaneously  but  separately.  This  is  accom¬ 
plished  by  adding  a  beam  splitter  and  a  second  analyzer/ 
photodiode  combination  to  the  reflected  light.  The  analyz¬ 
ers  are  adjusted  so  that  one  senses  the  component  parallel 
to  the  plane  of  incidence  and  the  other  detects  the  compo¬ 
nent  perpendicular  to  the  plane  of  incidence.  In  this  way, 
one  does  not  rotate  the  analyzer  to  detect  the  individual 
components. 
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ABSTRACT: 


The  magnetic  properties  of  molecular  beam  epitaxy  iron  films  grown  on  (001) 
GaAs  substrates  were  studied  using  a  SQUID  magnetometer.  In  uncapped  films, 
where  there  is  an  oxidized  iron  layer  on  top  of  the  film,  the  M-H  loops  shift 

from  being  symmetric  about  the  origin  -when  the  film  is  cooled  in  applied 

fields  to  temperatures  below  1G0K.  The  observed  behavior  is  attributed  to  a 
unidirectional  exchange  anisotropy  generated  by  the  antiferromagnetic  ordering 
of  the  oxidized  iron  surface .  By  comparing  the  magnitude  of  the  exchange 
coupling  observed  in  these  films  with  that  observed  in  the  Co-CoO  system,  the 

oxide  which  forms  on  the  iron  surface  is  most  likely  FeO.  In  addition,  we 

have  studied  the  magnetic  training  effect  or  the  effect  of  cycling  through  the 
hysteresis  loops  at  fixed  temperatures.  It  is  noted  that  the  temperature 
dependence  of  the  measured  exchange  coupling  is  different  from  that  which  is 
usually  reported.  The  reason  for  this  difference  is  unknown. 


PACS  number: 


75.30.Et,  75 . 50 . Bb ,  75.70. -i. 


Introduction 


Exchange  coupled  ferromagnetic  (F)  and 
antiferromagnetic  systems  demonstrate  a  shifted  hysteresis 
loop  when  the  systems  are  cooled  in  a  magnetic  field  to 
below  the  Neel  temperature  of  the  antiferromagnetic  (AF) 

layer. ^  Another  feature  ^f  this  phenomena  is  that  there  is 

2 

a  magnetic  training  effect,  i.e.,  at  a  fixed  temperature, 
the  magnitude  of  the  observed  effect  depends  upon  the 
magnetic  history  of  the  sample  at  that  temperature.  Although 
many  features  of  this  effect  are  not  well  understood,  we 
have  used  a  comparison  of  the  magnitude  of  the  shift 
observed  in  epitaxed  iron  films  to  that  in  Co-CoO  exchange 
coupling  to  identify  the  oxide  which  forms  on  the  epitaxed 
iron  film  surfaces  as  FeO.  In  addition,  the  temperature 
dependence  of  the  shift  and  the  coercivity,  and  its 
thickness  dependence  are  reported. 

Experimental  Details  and  Results 

The  iron  films  used  in  this  study  ware  grown  by 
molecular  beam  epitaxy  on  (001)  GaAs  substrates  at  200  C. 
After  film  growth,  the  normal  procedure  was  to  move  the  film 
to  the  introduction  chamber  after  the  film  had  cooled.  This 
chamhe-  was  then  filled  with  air  which  resulted  in  the 
formation  of  a  self -passivating  iron  oxide  layer,  removing 


As  a  control,  one  sample  was 


about  1.5  nm  thick  of  Fe. 
covered  with  amorphous  GaAs  in  the  growth  chamber  prior  to 
its  exposure  to  air,  prohibiting  the  iron  oxide  growth. 

The  magnetic  measurements  were  taken  following  a  number 
of  different  magnetic  field- temperature  paths.  For  most  of 
the  data,  a  magnetic  field  on  the  order  of  1000  Oe  was 
applied  in  the  plane  of  the  sample  along  a  magnetically  easy 
axes  with  the  sample  at  room  temperature,  cooled  to  a  lower 
temperature  and  a  hysteresis  loop  was  measured.  Next  the 
sample  was  warmed  to  a  higher  temperature  and  another 
hysteresis  loop  was  measured.  Other  data  were  accumulated 
by  always  warming  the  sample  to  room  temperature,  reapplying 
the  magnetic  field  and  then  cooling  to  another  temperature. 
In  addition,  for  measurements  made  when  the  samples  were 
cooled  in  zero  applied  field,  the  samples  were  always 
demagnetized  prior  to  cooling. 

The  main  experimental  observation  in  the  data  is  a 
shift  in  the  hysteresis  loop  for  the  iron  oxide  covered 
films  cooled  in  the  presence  of  an  applied  magnetic 
field  as  illustrated  in  fig.l.  This  shift  results  from  a 
coupling  between  the  ferromagnetic  iron  film  and  the 
antiferromagnetic  iron  oxide  which  forms  on  the  free  surface 
of  the  iron  films.  The  reason  to  attribute  the 
antiferromagnetism  to  the  oxide  is  that  the  data  taken  on 
the  sample  which  was  prepared  with  the  GaAs  protective 
overcoat,  therefore  lacking  an  iron  oxide,  did  not  exhibit 


any  features  indicative  of  exchange  coupling  at  any 
temperature . 

For  all  the  field  cooled  samples ,  a  measurable  shift  in 
the  hysteresis  loops  occurred  starting  at  temperatures  on 
the  order  of  100K  or  below.  Tne  data  were  quantified  in  the 
following  manner,  the  shift  from  the  origin  for  the 
hysteresis  loops  was  defined  as  the  exchange  field,  H  and 

one -half  the  total  width  of  the  loops  was  defined  as  the 

coercive  field,  H  .  For  all  the  data  shown  here,  the 
c 

samples  were  cooled  in  a  magnetic  field  to  a  low 
temperature,  usually  on  the  order  of  5K  and  then  the 
hysteresis  loops  were  measured  warming  the  sample  and 
stabilizing  the  temperature  for  each  loop.  The  fact  that 
there  is  a  training  effect,  which  wrill  be  discussed  next, 
introduces  some  error  in  the  data  accumulated  in  this 
manner,  however  by  using  both  techniques  described  above,  it 
was  determined  that  the  error  is  minimal.  The  data  showing 
the  temperature  dependence  of  He  and  for  the  four  samples 

measured  which  exhibit  this  behavior  are  shown  in  figures  2 
and  3 . 

Tne  magnetic  training  effect  refers  to  the  gradual 
reduction  in  the  exchange  field  and  coercivity  with  the 
number  of  hysteresis  loop  cycles  at  a  fixed  temperature. 

4 

Behavior  of  this  type  has  been  previously  reported  with 

2 

attempts  at  a  theoretical  understanding  ,  however  at  the 


present  there  does  not  appear  to  be  a  wholly  satisfactory 
model  for  this  effect.  One  feature  of  the  present  work  is 
there  appears  to  be  a  strong  thickness  dependence  to  the 
magnitude  of  the  training  effect.  Only  two  of  the  four 
samples  were  studied  in  depth  for  the  magnetic  training 
effect  but  of  these  two  samples,  the  training  effect  was 
considerably  larger  in  the  thinner  sample.  As  an  example, 
at  5K,  the  reduction  in  after  one  hysteresis  cycle  was 

45%  for  a  sample  2.0  nm  thick,  and  only  15%  for  a  sample  6.0 
run  thick. 

Discussion  of  Results 

2 

Simple  considerations  show  that  the  loop  shift  and  the 
ferromagnetic -antiferromagnetic  exchange  constant  are 
related  by  the  following  expression, 

W  -^K-  (1> 

where  tp  is  the  thickness  of  the  ferromagnetic  layer,  is 

the  loop  shift,  M  is  the  magnetization  of  the  ferromagnetic 

layer,  n  is  the  areal  density  of  spins  at  the  interface,  and 
J,,  is  the  exchange  constant  between  the  F-AF  spins  per  unit 

is. 

interface  area.  From  (1),  the  loop  shift  is  inversely 
proportional  to  the  thickness  of  the  ferromagnetic  layer. 
This  is  consistent  with  the  data  which  are  shown  in  fig.  4 


which  is  a  plot  of  the  H  's  measured  at  5K  versus  the 

e 

inverse  of  the  thickness  for  the  four  films  studied.  Using 
expression  (1) ,  the  magnitude  of  the  exchange  coupling 
constant,  J,, ,  between  the  F  and  AF  spins  across  the 

interface  was  found  to  be  9.3  X  10  ^erg  at  a  temperature  of 
5K.  For  comparison,  in  the  Co-CoO  system  the  zero 

temperature  value  of  J  was  found  to  be  5.59  X  10  ^ergs.~* 

K 

It  is  interesting  to  speculate  on  the  magnitude  of  the 
coupling  energy  one  should  expect  to  observe.  Although  it 
is  difficult  to  determine  the  interface  exchange,  one  might 
expect  it  to  be  related  to  the  exchange  energy  of  the  metal 
and  the  oxide.  Both  Co  and  Fe  have  comparable  exchange 
energies  so  one  might  expect  the  above  energy  differences  to 
be  associated  with  the  oxide.  In  the  case  of  CoO,  the  Neel 
temperature  is  on  the  order  of  room  temperature. 

Postulating  that  FeO  is  the  oxide  which  forms  on  the  iron 
films,  it  has  a  Neel  temperature  on  the  order  of  200K.  The 
ratio  of  the  two  Neel  temperatures  would  indicate  that  the 
value  expected  for  the  Fe-FeO  system  should  be  on  the  order 

of  3.8  x  10  ^ergs  which  is  in  reasonable  agreement  with  the 
experimental  result  given  that  the  above  argument  does  not 
include  any  differences  in  the  coupling  of  the  metals  to  the 
oxides . 

Tnere  has  been  a  report^  that  FeO  forms  on  iron  films 
in  an  oxygen  environment.  This  same  study  also  stated  that 


for  iron  oxides  grown  in  an  ambient  environment,  the  oxide 

which  forms  is  probablv  either  Fe.O.  or  7-Fe.O..  Vie  cannot 
‘  '  J  4  l  S 

positively  rule  out  the  possibility  of  either  one  of  these 

last  two  oxides  being  present,  however,  the  effect  which  is 

observed  must  have  its  origin  in  an  antiferromagnetic 

material.  It  is  certainly  possible  that  all  three  species 

are  present  which  might  also  explain  other  features  of  the 

data  which  are  discussed  below. 

2 

It  is  usual  for  the  to  be  linear  in  temperature. 

However,  it  is  clear  from  fig.  1  that  this  is  not  the  case 
for  the  samples  reported  here.  This  nonlinear  decrease  in 
Hg  might  come  from  the  temperature  dependence  of  the 

anisotropy  in  the  antiferromagnetic  iron  oxide.  In  the 

above,  we  have  assumed  that  the  antiferromagnetic  lattice  is 

locked  or  is  not  free  to  rotate  with  the  iron.  If  this  is 

not  the  case  then  the  H  which  would  be  measured  would 

e 

reflect  the  antif erromagnet  anisotropy  energy  and  not  the 
exchange  coupling.  As  an  example  of  this,  in  the  Ni-NiO 
system,  there  is  no  effect  of  the  antiferromagnetic  ordering 
of  the  NiO  on  the  hysteresis  loops  because  of  the  small 

value  of  the  anisotropy  energy  of  NiO.~  This  might  be  a 
possible  explanation  for  the  temperature  dependence  of 

and  ia  also  a  possible  explanation  why  the  J  measured  is 

K 

low  compared  to  the  Co-CoO  system.  On  the  negative  side  of 
this  argument,  if  this  is  the  case  then  one  might  expect  the 


training  effect  to  be  larger  than  normal.  This  possibility 
was  not  fullv  investigated,  however,  at  the  lowest 
temperatures  the  training  effect  was  comparable  to  that 

,  .  4 

ooserved  m  o~ner  systems. 

Summary' 

A  study  of  the  effects  of  exchange  coupling  in 
ferromagnetic  iron-antiferromagnetic  iron  oxide  has  been 
made.  As  expected,  the  magnitude  of  the  observed  shift  in 
the  hysteresis  loops  is  proportional  to  1/t  of  the  iron 
films.  By  comparing  the  magnitude  of  the  interfacial  energy 
in  the  Fe-Fe  oxide  films  to  that  observed  previously  in  Co- 
CoO  interfaces,  the  Fe  oxide  has  been  identified  as  FeO. 

The  temperature  dependence  of  the  exchange  field  is 
different  from  that  usually  observed.  This  may  be  due  to 
the  fact  that  the  FeO  has  a  Neel  temperature  of  200K  and 
therefore  over  the  temperature  range  investigated,  the 
anisotropy  energy  of  the  FeO  may  be  changing. 
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Figure  Captions: 


1.  The  measured  shifted  hysteresis  loop  at  51'  for  a  6 Or.  thick  file,  alter 
cooling  in  a  1000  Oe  applied  field. 

2.  The  measured  dependence  of  the  loop  shifts  or  H  on  temperature.  The 

symbols  plus  ~f-  ,  black  S  ,  cross  X  >  ant^  square  Q  correspond  to  sample 
thickneses  of  20A,  60A,  60A,  and  110A  respectively.  The  lines  are 
drawn  as  guides  to  the  eye. 

3 .  The  measured  coercivities  of  three  of  the  samples  as  a  function  of 
temperature.  The  sample  thicknesses  of  60A,  60A,  and  110A  correspond 
to  the  symbols  block  ■  ,  cross  X  .  and  square  □  respectively .  The 
lines  drawn  are  guides  to  the  eye. 

4.  The  magnitude  of  the  5K  loop  shifts  versus  the  inverse  of  the 
thickness . 
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Abstract 

For  many  materials  that  can  be  magnetized,  part  of  the  magnetization 
process  may  be  attributed  to  a  rotation  of  the  magnetization  vector.  In 
this  context,  a  combination  of  the  longitudinal  and  transverse  magneto¬ 
optical  Kerr  effects  are  used  to  detect  two  orthogonal  magnetization 
components  in  single -crystal  Fe/GaAs  (110)  thin  films.  Hysteresis  curves 
obtained  by  this  magneto-optical  technique  are  presented  for  fields  along 

the  in-plane  [001],  [liO],  and  [ill]  crystal  directions.  For  those  curves 
that  show  signs  of  rotation,  these  data  are  simulated  using  a  coherent 
rotation  mechanism  for  the  magnetization  process  and  Fresnel  reflection 
coefficients  for  the  two  Kerr  effects.  From  the  experimental  data,  it  is 

found  that  the  [111]  curves  have  shapes  that  are  indicative  of  a 

rotational  process.  On  the  other  hand,  both  the  [001]  and  [110]  have 

magnetization  curves  that  do  not  follow  a  simple  rotation.  From  the 

coherent  rotation  model,  there  is  qualitative  agreement  between  the 

modeled  and  experimental  data  for  the  [111], 


PACS  numbers:  78.20.Ls,  75.60.Ej 


Introduction 


For  ferromagnetic  materials  that  magnetize  partly  by  a 
rotation  process,  the  best  techniques  for  confirming  this 
rotation  utilize  the  vector  nature  of  the  magnetization. 
Within  this  context,  the  in-plane  magneto -optical  Kerr 
effects  can  be  utilized  to  detect  two  orthogonal 

1 

magnetization  components .  This  technique  is  used  to 
investigate  the  magnetization  processes  in  Fe/GaAs  (110) 
thin  films.  These  systems  have  been  extensively  studied  by 

2  3 

ferromagnetic  resonance,  ’  vibrating  sample  magnetometry 

*  s 

(VSM) ,  and  by  magnetoresistance  measurements.  The 

advantage  of  the  (110)  Fe  films  is  that  three  high  symmetry 

directions,  the  [001],  [110],  and  the  [111],  are  present  in 
the  plane  of  the  film. 

The  motivation  of  this  work  is  to  use  the  magneto¬ 
optical  technique  to  explore  the  rotational  processes  for 

applied  fields  along  the  [llO]  and  [ill].  Previous  VSM  data 
indicated  that  for  these  two  orientations,  the  shapes  of  the 

magnetization  curves  were  suggestive  of  a  rotational 
3 

process.  With  the  magneto-optical  technique,  it  will  be 
possible  to  track  the  direction  of  the  magnetization  vector 
during  the  rotation.  Thus,  for  those  directions  that  show 
definite  signs  of  rotation,  the  reversal  process  can  be 


simulated  with  a  coherent  rotation  model  for  the 
magnetization  process  and  Fresnel  coefficients  for  the  dual 
Kerr  effects.  In  this  way,  possible  modes  for  the 
rotational  process  can  be  eliminated  by  a  comparison  of  the 
modeled  and  experimental  data. 

Experimental 

The  films  used  in  this  study  were  prepared  with 

current  molecular  beam  epitaxy  (MBE)  techniques.  The  GaAs 

o 

(110)  substrate,  misoriented  by  1  to  the  [111],  was  etched 
in  a  standard  6:1:1  solution  of  H2S04 :H20:H202 .  The 

resulting  surface  oxides  of  the  substrate  were  removed  by 
o  o 

heating  to  620  C  in  vacuum.  A  700  A  GaAs  buffer  layer  was 
o 

then  grown  at  540  C  at  a  rate  of  5sec/laver  using  a  large  As 

.6 

overpressure  of  4x10  Torr.  This  large  As  overpressure  was 
used  to  prevent  Ga  puddling.  After  allowing  the  pressure  to 

.9  o 

fall  to  3x10  Torr,  100  to  150A  of  Fe  were  deposited  at  a 

o 

substrate  temperature  of  200  C.  Subsequent  RHEED 
measurements  revealed  good  epitaxial  growth  of  the  Fe  but  a 
rather  rough  surface. 

The  magneto-optical  measurements  were  made  at  room 

o 

temperature  using  a  HeNe  laser  at  A-6328A.  The  light 
incident  upon  the  film  was  polarized  in  the  plane  of 


incidence.  The  sample  was  mounted  on  a  rotatable  stage 
between  the  pole  faces  of  a  magnet  in  such  a  manner  that  the 
applied  field  was  perpendicular  to  the  plane  of  incidence 
formed  by  the  incident  and  reflected  light. 

The  detection  of  two  magnetization  components , 

individually  or  simultaneously,  depends  upon  the  analyzer 
l 

orientation.  basically,  for  analyzer  angles  near 

extinction,  the  component  parallel  to  the  plane  of  incidence 

o 

is  dominant  while  for  angles  90  from  extinction,  the 
component  perpendicular  to  the  plane  of  incidence  is 
important.  At  angles  intermediate  between  these  two 
extremes,  both  components  are  simultaneously  detected. 

Results  and  Discussion 

Depicted  in  Fig.  1  are  the  magnetization  loops  for  an 

applied  field  along  the  [001],  [llO],  and  [111]  with  an 

o 

analyzer  angle  of  45  from  extinction.  The  square 
magnetization  curve  for  the  [001]  is  indicative  of  the 
absence  of  any  rotation  beyond  the  transition  at  250  Oe . 

This  is  implied  by  the  constant  intensity  for  applied  fields 

in  excess  of  250  Oe .  In  contrast,  the  [ 110 ]  and  [111]  have 
shapes  that  can  be  interpreted  as  a  result  of  a  rotation. 

To  cross  check  this,  one  can  monitor  the  component 


perpendicular  to  the  applied  field.  Fig.  2  portrays  the 

magnetization  curves  for  the  same  crystal  directions  but 

o 

with  the  analyzer  now  rotated  to  within  5  of  extinction. 
Recall  that  for  this  orientation  of  the  analyzer,  the 
magnetization  component  parallel  to  the  plane  of  incidence 
is  sensed  to  a  greater  degree  than  that  in  Fig.  1. 
Considering  these  data,  in  addition  to  Fig.  1,  it  appears 
that  a  rotation  process  occurs  only  when  the  applied  field 

is  along  the  [111]  and  not  in  the  other  directions.  These 
conclusions  are  deduced  by  the  following  considerations. 
Comparing  figures  la  and  2a,  and  figures  lb  and  2b,  only  the 
saturation  to  saturation  intensity  dilrerence  is  altered, 
and  this  is  due  to  the  simple  crossing  of  the  polarizer  and 
analyzer.  A  comparison  of  figures  lc  and  2c  indicate  a 
qualitative  difference  in  the  shape  of  these  hysteresis 
curves  suggesting  the  presence  of  a  component  perpendicular 
to  the  applied  field. 

For  the  [ill]  curve,  assuming  the  magnetization  rotates 
toward  the  [001]  as  the  applied  field  is  reduced  from  900 
Oe ,  the  component  perpendicular  to  the  applied  field 
increases  as  the  applied  field  decreases,  as  is  suggested  by 
the  increasing  intensity  of  the  magnetization  curve.  At 
saturation,  either  at  positive  or  negative  900  Oe ,  the 
magnetization  lies  parallel  to  the  applied  field  resulting 


in  the  absence  of  a  component  perpendicular  to  the  applied 
field.  This  is  the  reason  why  the  intensity  is  so  similar 
at  the  two  saturated  states.  On  the  other  hand,  for  the 

[001]  and  [110],  only  the  small  change  due  to  the  switching 
of  the  component  parallel  to  the  applied  field  is  detected, 
no  analyzer  dependence  of  the  shape  of  the  curve  is  noted. 
This  suggests  that  for  these  two  directions,  beyond  the 
intial  transition,  a  rotation  of  the  magnetization  is  not 
present,  at  least  to  the  sensitivity  of  this  technique. 

To  model  the  magnetization  curves ,  two  ingredients  are 
required:  the  intensity  expression  that  relates  the 
photodiode  current  to  the  direction  of  the  magnetization  and 
the  polarizer/ analyzer  angles ,  and  the  magnetization 
process,  in  this  case  a  coherent  rotation  model. 

First,  the  photodiode  current,  which  is  proportional  to 
the  light  intensity  after  the  analyzer,  can  be  shown  to  have 
1 

the  form 


2  1  2  t 
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where  6  is  the  analyzer  angle  measured  relative  to  the 
plane  of  incidence,  m^  and  m^  are  the  components  of  the 


magnetization  perpendicular  and  parallel  to  the  plane  of 

incidence,  and  the  terms  of  the  form  r^  are  the  Fresnel 

P£ 

reflection  coeffients  for  the  transverse  or  longitudinal 
Kerr  effects , 

6 

Secondly,  previous  work  by  magnetoresistance  suggests 
that  these  systems  car.  be  modeled  using  a  ohenomenoiogical 
expansion  of  the  free  energy  density  of  the  form 

22  22  22  2  2  _ 

E  “  Kj  (QjQj  +  QjCtj  +  q2qs)  +  K^sin  6  +  2?rM^  -  M*H,  (2) 

where  in  this  expression  Kx  is  the  fourth  order 
magnetocrystalline  anisotropy  with  the  a^'s  being  the 

direction  cosines  of  M  in  a  cubic  system  such  as  Fe ,  is  a 

uniaxial  anisotropy  with  8  being  the  angle  that  M  makes  with 
the  [001];  this  term  could  be  due  to  strain  or  preferential 

2  7  2 

orientation  of  defects  in  the  Fe  film,  ’  2mM  is  the 

n 

demagnetization  energy  for  the  component  of  M  perpendicular 
to  the  film  plane,  and  the  last  term  is  the  energy  of  the 
magnetization  in  an  applied  field.  By  virtue  of  the  large 
demagnetizing  energy,  the  magnetization  is  effectively 
forced  to  lie  in  the  film  plane. 


For  a  nonzero  applied  field,  Eq .  2  is  solved 
numerically  for  the  stable  orientation  of  the  magnetization. 
The  magnetization  loop  is  the  result  of  tracking  the  local 
energy  minimum  as  a  function  of  H.  The  components  parallel 
and  perpendicular  to  the  plane  of  incidence  can  then  be 
determined  from  the  orientation  of  the  magnetization. 
Simulations  of  the  magnetization  processes  in  Fe/GaAs  (100) 

8 

films  by  this  method  resulted  in  agreement  with  the  data. 

o 

Depicted  in  Fig.  3  is  a  comparison  of  the  0^— 85  for 

the  simulated  and  experimental  data.  In  the  modeled 
magnetization  curve,  Fig.  3a,  K1/Ms~235  Oe ,  K^/M^-30  Oe ,  and 

o 

the  applied  field  is  misoriented  from  the  [111]  by  1  toward 
the  [001],  This  misorientation  was  used  to  induce  a 
rotation  in  the  assumed  direction.  These  anisotropy 

constants  that  fit  these  data  are  in  good  agreement  with 

6 

those  deduced  by  Riggs  et  al.  for  similar  films.  The 
simulated  curves  have  a  good  fit  to  the  data. 

For  the  film  prepared  here,  the  first  order  phase 
transition  was  not  observed  in  the  magnetization  data  for 

4  9 

the  [110],  ’  This  could  be  a  result  of  the  surface 
roughness  or  the  misorientation  of  the  (110)  GaAs  toward  a 

6  9 

[111].  For  the  anisotropy  constants  specified  by  others, 
the  simulations  here  predict  the  general  shape  seen  in  the 


VSM  data,  including  the  discontinuities  in  the 
magnetization  as  a  function  of  the  applied  field.  However, 
the  shape  of  the  simulated  curves  are  sensitive  to  the  value 
of  K  .  Thus,  if  the  surface  roughness  or  the  misorientation 

of  the  (110)  GaAs  can  induce  a  significant  change  in  the 

value  of  K  ,  the  absence  of  the  first  order  transition  could 
u 

be  accounted  for  in  this  film.  A  small  misorientation  of 

the  applied  field  from  the  [ 110 ]  does  not  appear  to  be  a 
plausible  explanation  since  the  discontinuity  is  still 

o  « 

present  even  for  misorientations  of  7  .  Another  feasible 
explanation,  given  the  large  corercivity  of  this  film  along 
the  [001]  and  the  equivalent  technical  saturation  field  for 

the  [110],  is  that  the  unpinning  energy  of  the  domains  is 
comparable  to  the  anisotropy  energy  resulting  in  p 

transition  to  the  [110]  without  ever  seeing  the  [111]  energy 
maximum. 

In  summary,  the  magneto- optical  technique  indicates 
that  the  magnetization  process  for  applied  fields  along  the 

[ill]  is  partly  due  a  rotation  process.  The  simulated 
curves  reproduce  the  analyzer  dependence  seen  in  the  data 

for  the  [111],  However,  the  [llO]  magnetization  curve  does 
not  show  signs  of  a  first  order  phase  transition  seen  in 


oj.mij.arly  prepared  films.  A  possible  explanation  could  be 
the  presence  of  a  rough  surface  or  the  misorientation  of  the 
substrate . 
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Figure  Captions 


Fig.  1  Magnetization  curves  for  an  applied  field  along  the  a.)  [001],  b .  ; 

[ilOj i  and  c. )  [ 111 j .  Tne  analyzer  is  oriented  45  from  the  plane  of 

incidence.  The  data  is  scaled  to  bring  the  saturation  states  of  I(-M  ) 

s 

and  I(+M  )  to  relative  values  of  +1  and  -1,  respectively . 

Fig.  2  Magnetization  curves  for  an  applied  field  along  the  a.)  [001],  b.) 

o 

[110],  and  c.)  [Ill]  with  the  analyzer  oriented  85  from  the  plane  of 
incidence.  The  curves  are  multiplied  by  a  factor  of  30  relative  to  those 
in  Fig.  1. 

Fig.  3  A  comparison  of  the  a.)  simulated  magnetization  curves  to  the  b.) 

o 

experimental  curves  for  the  [111]  for  an  analyzer  angle  of  85  from  the 

plane  of  incidence.  For  these  simulations,  K,/M  -235  Oe ,  K  /M  -30  Oe ,  and 

1  s  u  s 

0 

the  applied  field  is  oriented  1  from  the  [111]  toward  the  [001], 
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Abstract 

Systems  whose  dynamics  are  described  by  a  quasilogarithmic  or  stretched 
exponential  time  dependence  are  usually  fitted  by  models  which  use  disorder  to 
create  a  distribution  of  relaxation  times.  Here  we  describe  a  model  which 
decays  slowly  towards  equilibrium  but  does  not  require  disorder  to  provide  the 
slow  dynamic.  The  model  is  consists  of  a  spin  system  with  the  spins 
interacting  via  the  dipole-dipole  interaction.  The  model  is  able  to  replicate 
the  more  pronounced  features  observed  in  the  magnetization  decay  of  magnetic 
systems  and  high  temperature  superconductors. 

PACS  numbers:  05.50.+q,  61.42.+h,  64.70.Pf,  75. 10. Nr 
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2 


A  number  of  physical  systems  are  known  to  relax  more  sloviv  than  expected 
fiom  a  simple  Debye  relaxation  model.  These  systems  include  glasses, 

1  2 

polymers  ,  and  spin  glasses.  Usually  a  distribution  of  relaxation  times,  a 
manifestation  of  disorder,  is  argued  to  be  the  origin  of  the  slow  relaxation 

3 

of  some  property  of  the  system  .  Within  the  context  of  a  distribution  of 
relaxation  times,  the  relaxing  entities  are  noninteracting,  i.e.,  the 
relaxation  process  occurs  in  parallel  for  all  the  entities.  Using  this 
generic  model  to  analyze  physical  systems  merely  yields  a  distribution 
function  describing  the  disorder  which  may  or  may  not  be  identifiable  with 
some  physical  property  of  the  system.  Similarly,  computer  simulations  have 
also  been  plagued  by  their  inability  to  make  connections  to  real  physical 
systems.  Another  theoretical  approach  has  considered  serial  relaxation 
4- 

processes  rather  than  parallel  relaxation.  The  model  of  heirarchially 
constrained  dynamics  that  result  from  this  approach  do  contain  slow  relaxation 
but  are  again  phenomenological  in  nature. 

In  this  letter  we  present  a  surprisingly  simple  model  which  requires  no 
disorder  to  provide  the  slow  relaxation  observed  in  so  many  systems.  In  this 
model,  interactions  provide  for  the  slow  dynamic  which  are  manifest  by  either 
quasilogarithmic  or  stretched  exponential  behavior.  In  particular  the  model 
predicts  the  nonmonotonic  temperature  dependence  observed  in  the  decay  of  the 

remanent  magnetization  of  CoCr  films'*  and  high  temperature  superconductors.^’^ 
Although  not  contained  in  the  present  paper,  it  may  be  used  as  a  direct 
comparison  of  the  dynamics  generated  by  Monte  Carlo  calculations  on  spin 


systems . 
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The  svstem  consists  of  a  planar  array  of  spins  with  an  anisotropy  energy 
which  is  a  minimum  when  the  spins  are  perpendicular  to  the  plane.  The  model 
includes  the  dipole-dipole  interaction  between  the  spins  which  are  ir.  essence , 
the  causal  agent  of  the  slow  relaxation.  Tne  interactions  are  treated  m  the 
mean  field  limit  using  the  demagnetization  field  arrising  from  the  plane  of 
spins . 

Due  to  the  dipole-dipole  interaction  between  the  spins,  the  ground  state 
possesses  no  net  magnetization  perpendicular  to  the  plane  even  though  the 
anisotropy  energy  for  a  single  spin  prohibits  the  spins  from  lying  in  the 
plane;  the  ground  state  is  in  effect  an  antiferromagnet .  If  the  system  is 
initially  polarized  with  all  the  spins  oriented  perpendicular  to  the  plane 
there  will  be  a  demagnetization  field  which  drives  the  system  to  the  ground 
state  or  the  zero  net  magnetized  state.  Because  of  the  anisotropy  energy' 
inhibiting  the  relaxation  process  and  the  fact  that  the  drive  field  reduces  as 
the  system  relaxes,  there  is  a  slow  demagnetization  of  the  spins.  In  what 
follows,  we  will  outline  the  model,  exhibit  some  of  its  interesting  features, 
and  discuss  how  it  may  pertain  to  a  variety  of  other  systems . 

Tne  fundamental  unit  of  the  model  is  a  single  magnetic  entity  such  as  a 
grain  of  a  magnetic  material  which  has  a  volume  v,  and  a  saturation 
magnetitation  of  m  .  Tnis  grain,  or  spin  as  we  refer  to  it,  has  a  magnetic 

anisotropv  of  the  form  -K  cos  ( 6 )  where  is  positive  and  6  is  the  angle 

between  the  magnetization  vector  and  the  +z  axi s.  In  the  presence  of  an 
external  field  H,  in  the  -z  direction,  the  energy  of  the  magnetic  grain  is 
given  by 

E  -  -K  v  co s2(0)  +  K  Hv  eos(0). 
ti  s 


(1) 
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This  equation  is  the  starting  point  for  Neel's  nndel  rr  flue  grains  usee  to 
model  the  temporal  evolution  the  magnetisation  of  the  solidified  magma  ir.  the 
earth's  crust. 

Of  particular  interest  to  the  present  discussion  is  to  consider  the 
relaxation  of  a  number,  N,  of  the  above  particles  which  lie  in  the  x-y  plane 
in  a  sheet  one  layer  thick.  The  N  particles  or  spins,  as  they  will  be  refered 
to  usually,  are  prepared  in  a  state  with  the  magnetization  saturated  in  the  +z 
direction.  The  process  of  reversing  the  magnetization  of  a  particle  requires 
the  rotation  of  the  magnetization  past  a  energy  barrier  which  is  due  to  the 
anisotropy  energy.  In  the  presence  of  a  reversed  magnetic  field  (-z 
direction) ,  one  may  determine  the  value  of  6  for  the  maximum  value  of  the 
energy  by  differentiating  eq . 1  with  respect  to  6  and  setting  the  derivative 
equal  to  zero.  Following  this  procedure  one  finds  the  energy  maximum  ED , 

separating  the  -t-z  and  -z  magnetic  states  to  be  given  by 

E  -<m  2H2v)/(4K  .).  (2) 

D  S  U 

In  order  to  determine  the  dynamics  it  is  the  energy  differences  between  the  up 
and  down  states  relative  to  the  energy  maximum  or  barrier  which  are  relevant. 
These  energies  are  given  by 

E+(.)-tv/(4Ku)](msH-(+)2Ku)2 


(3) 
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where  the  ^  refers  to  the  energy  difference  between  the  up  (down)  state 

and  the  energy  barrier.  Using  these  energies,  the  rate  equation  is  then  giver, 
by 


,  ..  -0E+  -0E- 

dn  /dt-  -w  n  e  +w  n  e 
+'  +  +  -  - 


(4) 


where  n+  is  N+  ^ .  is  the  number  of  antiparallel  (parallel)  divided  by  N, 

P  is  l/k^T  and  w+^  is  the  attempt  frequency  from  the  up  (down)  state, 
antiparallel  to  the  applied  field,  to  the  down  (up)  state.  Expressions  for 

q 

the  w's  which  depend  upon  the  state  of  the  system  have  been  derived  by  Brown 

and  by  Neel^  but  our  calculations  indicate  that  very  little  inaccuracy  is 
induced  by  considering  the  w's  to  be  equal  and  independent  of  the  applied 
magnetic  field. 

If  the  decay  of  the  magnetic  system  described  by  equation  4.  is  monitored 
at  a  fixed  time  eg.  30  sec,  as  a  function  of  the  applied  magnetic  field,  one 
will  find  a  maximum  in  the  decay  rate  when  the  applied  field  is  equal  to 
2Ku/ms ,  the  coercive  field  for  this  model  system.  Physically  this  can  be 

understood  by  realizing  that  this  field  corresponds  to  the  disappearance  of 
the  energy  barrier  separating  the  up  and  down  states.  For  fields  less  than 
this  field,  the  energy  barrier  is  too  large  to  allow  the  magnetic  moments  to 
flip  and  the  rate  dn+/dt  is  small.  When  the  applied  field  is  larger,  most  of 

the  magnetization  has  flipped  so  dn+/dt  is  small  at  the  prescribed  time.  It 

is  only  when  the  applied  field  is  equal  to  the  coercive  field  that  the  decay 
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rate  from  +  to  -  is  large . 
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Vie  now  consider  the  modification  to  this  model  which  provides  for  the 

slow  dynamics;  allow  the  spins  to  interact  via  the  dipole-dipole  interaction. 

As  before,  consider  the  N  grains  of  the  system  to  be  confined  to  a  plane  with 

the  easy  direction  of  the  magnetization  of  the  grains  to  be  perpendicular  to 

the  plane.  In  the  mean  field  situation  the  dipole -dipole  interaction  is  given 

by  the  demagnetization  field  of  the  plane  of  spins  interacting  with  the 

grains.  The  demagnetization  field  is  given  by  4xM  where  M,  the  magnetization 

of  the  sample  is  given  by  m  (n  -n  ) .  By  defining  a  fin  such  that  n  -1/2  +  fin/2 

s  +  + 

and  n  —1/2+  fin/2  then  the  demagnetization  field  is  just  47rms5n. 

If  the  planar  sample  of  magnetic  grains  is  first  polarized  with  all  the 
magnetization  in  the  +z  direction,  then  the  demagnetization  field  is  directed 
opposite  the  magnetization  and  will  tend  to  drive  the  system  into  a 
demagnetized  state.  But  as  the  magnetization  of  the  system  decays,  the 
driving  field  or  the  demagnetizing  field  will  also  become  smaller  thereby 
reducing  or  slowing  down  the  dynamics  of  the  system.  This  is  the  essence  of 
the  origin  of  the  quasilogarithmic  decay  of  the  magnetization  of  the  system. 

If  the  applied  field  in  eq.3  is  replaced  by  the  demagnetizing  field,  the 
above  definitions  are  used,  and  after  a  little  algebra,  then  the  rate  of  the 
decay  of  fin  is  given  by 


dfin/dt  -  -w+(l/2+6n/2)exp{(-^vKu)[(2»rMs2fin/Ku)-l]2}  + 
w_ (l/2-5n/2)exp{ (-0vK  ) [ (2wMs25n/Ku)+132} .  (5) 

Again  for  this  situation  we  should  note  that  one  may  derive  relationships  for 
the  w's  but  in  solutions  of  eq.5  little  difference  is  exhibited  if  the  w's  are 
equal  and  constant.  Equation  5  describes  the  relaxation  of  a  system  of 
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particles  where  the  drive  force  for  the  relaxation  depends  upon  the 
instantaneous  state  of  the  system.  As  stated  earlier,  this  expression  is  in 
the  mean  field  limit  and  does  not  consider  the  actual  discreteness  of  the 
relaxing  grains  nor  their  particular  geometry.  This  expression  is  amenable  to 
numerical  techniques  for  its  solution  with  different  values  of  the  parameters 
in  the  expression.  Figure  1  shows  an  example  of  eq.5  solved  using  a  Runge- 
12 

Kutta  algorithm  and  indicates  that  the  relaxation  is  quasilogrithmic  over  7 
decades  in  time!  The  recovery  of  a  simple  Debye  relaxation  of  this  model  is 
accomplished  when  the  argument  of  the  exponentials  are  less  than  unity  which 
occurs  when  the  magnetic  entities  are  in  the  superparamagnetic  limit.  A 
useful  or  effective  approximation  to  this  limit  is  >  vK^;  this  is  the  usual 

condition  for  the  superparamagnetic  limit. 

Two  systems  which  would  appear  to  be  ideally  modeled  by  the  above  are 
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saturated  films  of  CoCr  which  have  a  perpendicular  anisotropy  and 
superconductors  which  have  been  cooled  in  a  magnetic  field  trapping  flux.  In 
both  cases  the  systems  are  out  of  equilibrium  possessing  a  net  magnetization 
thereby  providing  a  decaying  driving  force.  Also,  their  dynamics  may  be 
described  by  incremental  changes  in  the  magnetic  state,  i.e. ,  for  the  CoCr 
films  the  magnetization  in  a  single  grain  flips  its  magnetization  direction, 
and  in  the  superconductor  a  flux  quantum  leaks  out  of  the  sample. 

In  the  case  of  CoCr  films,  the  morphology  of  the  alloy  consists  of  small 

14 

columnar  grains  of  a  magnetic  alloy  imbedded  in  a  nonmagnetic  matrix.  The 
column  axes  of  the  grains  are  perpendicular  to  the  plane  of  the  film  and 
define  the  magnetically  easy  direction  for  the  magnetic  moment  of  the  grains. 
In  the  superconductor  case,  it  is  again  the  remanent  magnetization  which  is 

describable  by  the  above  model. ^  When  the  magnetic  field  is  removed,  the 
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superconductor  will  generally  maintain  some  density  of  vortices  pinned  by 
detects  which  gives  a  situation  analogous  to  the  magnetic  grains  discussed 
above.  The  pinning  energies  play  the  role  of  the  anisotropy  energy  in  the 
magnetic  film. 

We  should  point  out  one  rather  strong  similarity  of  the  two  systems  and 
their  remanence  decay.  All  the  simple  models  which  give  a  slow  decay,  either 
logarithmic  or  quasilogrithic ,  also  predict  a  monotonic  increase  in  the  decay 
rate  with  increasing  temperature.  In  CoCr  films  and  in  high  temperature 
superconducting  samples  the  experimental  evidence  is  contrary  to  this 
prediction;  in  both  cases  the  decay  rate  passes  through  a  maximum  value  with 
increasing  temperature . 

The  model  presented  also  contains  such  behavior.  At  low  temperatures  the 
remanent  magnetization  is  large  but  because  there  is  little  thermal  energy  the 
decay  slope  is  small.  At  very  high  temperatures  the  thermal  energy  is 
sufficient  for  the  system  to  relax  rapidly  but  in  the  time  window  of  the 
measurement  most  of  the  magnetization  has  decayed  and  therefore  the  decay 
slope  is  reduced.  This  is  very  similar  to  the  rationale  of  why  there  is  a 
peak  in  the  decay  rate  at  the  coercive  field  in  the  Neel  model  described 
previously.  This  nonmonotonic  behavior  is  shown  in  the  inset  of  fig.l. 

One  should  note  that  the  argument  of  the  exponential  in  Eq .  5  is 
quadratic  in  the  decaying  variable.  This  arises  from  the  way  the  applied 
field  determines  the  6  for  the  maximum  energy.  In  a  different  system  which 
does  not  behave  in  this  manner,  the  argument  of  the  exponential  may  be  linear. 
In  this  case  eq  5,  again  gives  a  function  which  decays  quasilogarithmically 
for  long  times.  The  linear  form  may  be  more  relevant  for  systems  where  the 
free  energy  provides  the  driving  force  i.e.,  the  argument  of  the  exponential 
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is  of  the  fcrir.  dF/kpT  where  dF  is  the  difference  ir,  free  energy  of  the 

instantaneous  configuration  from  equilibrium. 

For  application  to  other  systems,  such  as  glasses,  one  mav  object  to  the 
use  of  a  mean  field  approximation.  In  the  case  of  structural  glasses  the 

3 

strain  fields  are  long  range  falling  ofr  as  1/r  just  as  the  dipole-dipole 
interaction;  therefore  a  mean  field  approximation  should  be  applicable.  In 
the  case  of  spin  glasses  the  dipole-dipole  interaction  is  expected  to  play  an 
important  role  just  as  in  the  magnetic  model  presented  here. 

A  final  point  to  make  is  the  connection  between  the  present  model  and 

U 

that  of  the  theoretical  model  which  generates  a  stretched  exponential  via  a 
hierarchy  of  constraints.  In  that  paper,  with  certain  approximations,  is 
derived  a  stretched  exponential  relaxation  of  the  form 

M(t)-M0exp[-(t/r/],  O<0<1.  (6) 

As  argued,  this  form  of  relaxation  is  often  observed  in  strongly  interacting 
systems.  Within  the  present  model,  we  have  been  unable  to  fit  this  form  of 
the  relaxation  over  the  entire  time  domain  of  the  relaxation  process,  but  do 
find  our  model  can  be  approximated  by  the  stretched  exponential  form  over 
several  decades  in  time,  as  shown  in  figure  2.  This  is  the  same  data  as 
exhibited  in  figure  1  but  plotted  in  a  magnetization  domain  where  the 
stretched  exponential  form  applies.  In  the  plot  as  shown,  a  stretched 
exponential  decay  should  be  a  straight  line.  The  question  as  to  whether  or 
not  the  stretched  exponential  form  should  be  applicable  for  all  times  m  a 
real  physical  system  remains  open,  but  in  at  least  spin  glasses  it  does  not 
appear  to  hold.  Consequently,  the  fact  that  our  model  can  only  be  fitted  by  a 
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stretched  exponential  over  a  limited  time  domain  is  not  a  serious  deficit  of 
the  model. 

Acknowledgements 

One  of  the  authors  (EDO)  would  like  to  thank  Professors  J.  Kakalios ,  L. 
Glazman,  G.  Mazenko ,  0.  Vails,  and  E.  Shklovskii  for  a  number  of  interesting 
discussions.  This  research  was  supported  in  part  by  AFOSR  grant  AFOSR-89- 
0248. 


References 


I. K.L.  Ngai,  Comments  Solid  State  P'nvs .  9,  127  (1579), 
and  9 ,  141  (1979)  . 

2.O.S.  Lutes  and  J.L.  Schmidt,  P'nvs.  Rev.  12  5 .  433 
(1961)  and  Y.  Ishikawa,  R.  Tourr.ier,  and  J.  PnilipDi,  J. 
Phys.  Chem.  Solids  21.  1727  (1965). 

3. J.J.  Prejean  and  J.  Souletie,  J.  de  Physiaue  41,  1335 
(1980). 

4. R.G.  Palme  D.L.  Stein,  E.  Abrahams,  and  P.W. 
Anderson,  Phys.  Rev.  Lett.  53.,  958  (1984). 

5. D.K.  Lottis,  E.  Dan  Dahlberg,  J.A.  Christner,  J.I. 

Lee,  R  L.  Peterson,  and  R.M.  White,  Jour,  de  Phys.  49 . 
C8-1989  (1989). 

6. K.  Tuomenen,  A.M.  Goldman,  and  M.L.  Mecartney,  Phys. 
Rev  BIZ,  548  (1988). 

7. C.W.  Hagen  and  R.  Griessen.  Phys.  Rev.  Lett.  62,  2857 
(1969). 

8. L.  Neel,  Rev.  Mod.  Phys.  25 .  293  (1953)  and  Chikazumi 
"Physics  of  Magnetism"  (Wile}',  New  York,  1972)  Chapter 
15. 

9.  William  Fuller  Brown,  Jr.,  Jour.  Appl.  Phys.  3(),  130S 
(1959). 

10.  L.  Neel,  Compt.  Rend.  228 .  99  (1949). 

II.  One  should  note  that  the  data  are  taken  in  a  fixed 
time  w'indow.  Assuming  that  the  relaxation  of  the  system 
is  not  to  fast  nor  to  slow  allows  the  peak  to  be  at  the 
coercive  field.  Although  it  would  appear  that  the 
window  restriction  would  invalidate  the  analysis,  in 
practice,  as  long  as  one  is  in  a  temperature  limit  where 
a  slow  dynamic  is  observed,  it  is  a  reasonable 
assumption. 

12.  W.H.  Press,  B.P.  Flannery,  S.A.  Teukolsky,  and  W.T. 
Vetterling,  "Numerical  Recipes,"  (Cambridge  University 
Press,  New  York,  1986)  pp.  550-554. 

13.  D.K.  Lottis,  R.M.  White,  and  E.  Dan  Dahlberg, 
submitted  for  consideration  of  publication. 

14.  M.R.  Kahn,  J.I.  Lee,  D.J.  Seagle ,  and  N.C. 

Femelius,  Jour.  Appl.  Phys.  63,  833  (1988). 


12 


Figure  1. 

An  example  of  a  slow  decay  of  the  model  system 
given  by  equation  £ .  For  this  decay  the  relevant 

parameters  are  K  -5xl0''ergs ,  M^-200  emu/cm-' ,  T-  300K, 
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v-2xl0  ,  and  v  — v  -2x10  sec  * .  Note  that  the  decav  is 

-r 

measured  over  nine  decades  in  time.  The  insert  in  the 
figure  is  the  temperature  dependence  of  the  logarithmic 
slope  of  the  decay  taken  at  a  time  of  100  sec..  The 
nonmonotonic  behavior  of  the  slope  is  observed  in 
magnetic  and  superconducting  systems  (see  refs.  5,6  and 
7). 

Figure  2. 

This  is  the  data  from  figure  1  plotted  to  show  its 
compatibility  with  a  stretched  exponential  relaxation 
model.  As  shown,  the  stretched  exponential  is  well 
described  by  this  system  for  over  4  decades  in  time. 
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